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Dust formation history of galaxies: a critical role of metallicity*for the dust mass 
growth by accreting materials in the interstellar medium 
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This paper investigate what is the main driver of the dust mass growth in the interstellar medium (ISM) by using 
a chemical evolution model of galaxy with metals (elements heavier than helium) in dust phase in addition to the 
total amount of metals. We consider asymptotic giant branch (AGB) stars, type II supernovae (SNe II) and the dust 
mass growth in the ISM as the sources of dust, and SN shocks as the destruction mechanism of dust. Further, to 
describe the dust evolution precisely, our model takes into account the age and metallicity (the ratio of metal mass 
to ISM mass) dependence of the sources of dust. We particularly focused on the dust mass growth, and found that 
the dust mass growth in the ISM is regulated by the metallicity. To quantify this aspect, we introduce a "critical 
metallicity", which is a metallicity at which the contribution of stars (AGB stars and SNe II) equals that of the dust 
mass growth in the ISM. If the star formation timescale is shorter, the value of the critical metallicity is higher, 
but the galactic age at which the metallicity reaches the critical metallicity is shorter. From observations, it was 
expected that the dust mass growth was the dominant source of dust in the Milky Way and dusty QSOs at high 
redshifts. By introducing the critical metallicity, it is clearly shown that the dust mass growth is the main source of 
dust in such galaxies with various star formation timescales and ages. The dust mass growth in the ISM is regulated 
by metallicity, and we stress that the critical metallicity works as an indicator to judge whether the grain growth in 
the ISM is dominant source of dust in a galaxy, especially because of the strong and nonlinear dependence on the 
metallicity. 
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1. Introduction 

Stellar light, in particular at shorter wavelengths, is ab- 
sorbed by dust and re-emitted as the far-infrared thermal 
emission from the dust (e.g., Witt & Gordon, 2000, and ref- 
erences therein). Therefore dust affects the spectral energy 
distributions of galaxies (e.g., Takagi et al., 1999; Granato 
et al., 2000; Noll et al., 2009; Popescu et al., 2011). The 
existence of dust in galaxies also affects the star formation 
activity. Dust grains increase the molecular formation rate 
by two orders of magnitude compared to the case without 
dust (e.g., Hollenbach & McKee, 1979), and the interstellar 
medium (ISM) is cooled efficiently by molecules and dust. 
Consequently, star formation is activated drastically by dust. 
Hence, dust is one of the most important factors for the evo- 
lution of galaxies (e.g., Hirashita & Ferrara, 2002; Yama- 
sawa et al., 201 1). 

The amount of dust in galaxies is one of the crucial factors 
to interpret the observational information of galaxies, since 
dust exists ubiquitously and the radiation from stars is al- 
ways affected by dust attenuation. However, in spite of its 
importance, the evolution of dust amount has not been com- 
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pletely established yet. There are some key factors to un- 
derstand dust evolution of galaxies. One of the keys is the 
ratio of metal (elements heavier than helium) mass to ISM 
mass, which is called "metallicity". Since dust grains con- 
sist of metals, it is natural to think that the evolution of dust 
is closely related to metallicity. In general, galaxies are be- 
lieved to evolve from the state with a very low metallicity 
and very small amount of dust to higher amounts of metal 
and dust. Hence, it is a mandatory to model the formation 
and evolution of dust grains in galaxies along with the evo- 
lution of metallicity (e.g., Dwek & Scalo, 1980; Hirashita, 
1999a,b; Inoue, 2003; Yamasawa et al., 201 1). 

Dust grains are formed by the condensation of metals. A 
significant part of the metals released by stellar mass loss 
during stellar evolution or supernovae (SNe) at the end of 
the life of stars condense into dust grains. Dust grains are 
not only supplied from stars but also destroyed by SNe blast 
waves (e.g., Jones et al., 1994; Jones, Tielens & Hollenbach, 
1996; Nozawa et al., 2003; Zhukovska, Gail & Trieloff, 
2008). In addition, we should consider the dust mass growth 
in the ISM by the accretion of atoms and molecules of re- 
fractory elements onto grains (e.g., Dwek, 1998; Liffman & 
Clayton, 1989; Draine, 2009; Jones & Nuth, 201 1). 

What kind of dust formation processes are dominant at 
each stage of galaxy evolution? It is a very important ques- 
tion for understanding the evolution history of the ISM and 
star formation in galaxies. However, since dust evolution 
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depends strongly on the age and metallicity of a galaxy, it 
is not easy to answer this question. Up to now, dust evolu- 
tion has been studied with various models. For example, in 
young galaxies, SNe have been considered as the source of 
dust because they are the final stage of massive stars whose 
lifetime is short, and asymptotic giant branch (AGB) stars 
have been neglected because of their longer lifetime. How- 
ever, Valiante et al. (2009) showed that the AGB stars also 
contribute to the dust production in young galaxies and can- 
not be neglected even on a short timescale of ~ 500 Myr. 
A more elaborate survey of the parameter space for the dust 
formation by SNe and AGB stars has been done by Gall et 
al. (201 la). They showed that the contribution of AGB stars 
exceeds that of SNe II at several 100 Myr if the ratio between 
metal and dust mass produced by SNe II is less ~ 0.01 and 
mass-heavy IMF with mass range 1-100 M©. 

As for the dust mass growth, the ISM is considered to be 
the main source of dust in various galaxies. For example, the 
present dust amount observed in the Milky Way cannot be 
explained if the source of dust would have been only stars, 
suggesting that we must consider the dust mass growth in the 
ISM in evolved galaxies (e.g., Dwek, 1998; Liffman & Clay- 
ton, 1989; Draine, 2009; Jones & Nuth, 2011). Recently, 
dusty quasars (total dust mass > 10 8 M©) have been discov- 
ered at high redshifts (e.g., Beelen et al., 2006; Wang et al., 
2008), and theoretical studies on dust sources at high red- 
shifts are currently carried out actively (e.g., Michalowski et 
al., 2010b; Gall et al., 2011a,b; Pipino et al., 2011; Valiante 
et al., 2011). They showed that it is hard to explain the to- 
tal dust amount in these QSOs only with stellar contribu- 
tions, and then they discussed the importance of the dust 
mass growth in the ISM. Next question is what controls the 
point where the dust mass growth in the ISM becomes effi- 
cient to the total dust mass in galaxies. Therefore, although 
each physical process has been already extensively discussed 
in preceding studies, there emerges a crucial question: what 
kind of dust production process is dominant at each stage 
of galaxy evolution? Particularly, when the dust dust mass 
growth becomes dominant as a source of dust mass? 

The central aim of this work is to address this question. In 
this paper, we investigate what is the main driver of the dust 
mass growth in the ISM. Since all sources of dust production 
are tightly related to each other on dust evolution, it is crucial 
to treat these processes in a unified framework to understand 
the evolution of dust in galaxies. Here, we adopt the model 
based on chemical evolution model in a same manner as Hi- 
rashita (1999b); Calura, Pipino & Matteucci (2008); Inoue 
(201 1). This is because their models consider main dust pro- 
duction/destruction processes that affect the dust evolution 
of galaxies, and it is easy to compare our results to the previ- 
ous ones. From this work, we find that the dust mass growth 
in the ISM is regulated by metallicity. We refer to this metal- 
licity as the critical metallicity, whose details are described 
in Sect. 3.2. Although the dust mass growth can occur at any 
time of a galaxy age, but we stress this point because there 
is a moment at which the dust mass growth overwhelms the 
contribution from other sources of dust. 

This paper is organized as follows. In Sect. 2, we describe 
the model developed for this work. In Sect. 3, we show 
and discuss the basic results obtained by our model. The 



main topic of this paper, critical metallicity, is introduced and 
extensively examined in Sect. 3.2. Section 4 is devoted to 
the conclusions. The solar metallicity is set to be Z Q = 0.02 
(Anders & Grevesse, 1989) throughout this paper. 

2. Dust evolution model of galaxies 

In this section, we show the simple chemical evolution 
model with dust which examine what determines the point 
where the dust mass growth in the ISM becomes effective. 
The dust evolution model is built with the same manner as in 
Hirashita (1999b) and Inoue (201 1). 
2.1 Equations of galaxy evolution 

In this subsection, we describe the equations of mass evo- 
lution of stars and the ISM which contains metal and dust in 
galaxies. We treat a galaxy as one-zone because we are in- 
terested in global properties of galaxies. Also, we assume a 
closed-box model. Thus, the total baryon mass M tot (the 
sum of the stellar mass and the ISM mass) is a constant. 
However, since M to t is just a scale factor in our model, this 
value does not affect the physical properties of galaxies non- 
linearly. 

In this work, we do not consider the effects of inflow and 
outflow. However, they may not influence the properties 
of dust and metal enrichment in galaxies because of the 
following reasons: An inflow makes not only metallicity but 
also dust-to-gas mass ratio small, because usually an inflow 
is considered to be metal- and dust-poor. As for an outflow, 
it expels ISM components (gas, metal and dust) out of a 
galaxy. However, if all ISM components flow out together, 
the metallicity and dust-to-gas mass ratio do not change. 

Under these settings, equations of the time evolution of the 
total stellar mass M_, ISM mass Mism, metal mass Mz, and 
dust mass M d are (e.g., Lisenfeld & Ferrara, 1998; Hirashita, 
1999b) 



dM.(t) 
dt 

dM lSM (t) 
dt 
dM z (t) 

dt 
dM d (t) 



= SFR(i) - R(t), 

= -SFR(t) + R(t), 

= -Z(i)SFR(i) + R z {t) + Y z (t), 



(1) 
(2) 
(3) 



dl TSN T acc 



(4) 



where SFR is the star formation rate, Z(t) = Mz/Mism is 
the metallicity, V = M d /MisM is the dust-to-gas mass ratio, 
S = Md/Mz is the fraction of the metals in dust. r\ is the 
mass fraction of cold clouds where the accretion process oc- 
curs, tsn and T acc are the timescales of dust destruction and 
accretion, respectively. The definitions of these timescales 
are described later. Also, R(t) and Rz {t) are the total baryon 
mass returned by stars and the total metal mass once injected 
in stars and just returned in the ISM per unit time, respec- 
tively. Yz(t) and Y<i(t) are the total metal mass newly pro- 
duced and ejected by stars and the total dust mass ejected by 
stars per unit time, respectively. There is another notation 
which Yz includes Rz (e.g., Inoue, 2011). In this case, the 
value of Yz is different from that in this paper because of the 
different definition. 
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For the SFR, we adopt the Schmidt law (Schmidt, 1959); 
SFR oc Mjg M . Here, we adopt n = 1 for simplicity. Thus, 
the SFR is expressed as 



SFR(t) = 



M ISM (t) 



TSF 



(5) 



where tsf is the star formation timescale. 

Also, R(t), R z (t), Y z (t) and Y d (t) are written by 

/•100 M 

R(t) = / [m-w(m,Z(t-T m ))} 

J "icut(t) 

x</>(™)SFR(i - r m )dm, (6) 

"100 M 



/■1UU IVlQ 

(t) = / [m-w(m,Z(t-r m ))] 

J TO cu t(t) 



X(/)(m)SFR(t - r m )Z(t - r m )dm, 

(7) 

<-100 M 

it) = / m z (m, Z(t-T m )) 

J m cut (t) 

x0(m)SFR(t - r m )dm, (8) 

piao m 

y d (t) = / m d (m,Z(t - r TO )) 

m cut (t) 

x ^(m)SFR(i - r m )dm, (9) 

where </>(to) is the initial mass function (IMF), 
uj(m, Z(t)),mz(m, Z(t j) and md(m, Z(f)) are the 
remnant mass, metal mass newly produced and ejected, and 
dust mass produced and ejected by a star of initial mass m 
and metallicity Z, respectively. The lifetime of a star of 
initial mass m is expressed as r m , is taken from Raiteri, 
Villata & Navarro (1996). In this work, since its metallicity 
dependence is weak, we always adopt the lifetime of the 
case of the solar metallicity. The lower limit m cu t (t) is the 
mass of a star with the lifetime r m = t. As for IMF, we 
adopt Larson IMF (Larson, 1998) in the stellar mass range 
0.1-100 M , 



(m) oc m ( Q+1 ) exp ^~ 



m ) 



(10) 



Here we adopt a = 1.35 and m c h = 0.35 M©. Also, we 
normalize it as 



100 Mr 



m<f)(m)dm = 1. 



(11) 



0.1 M ( 



To calculate above four equations [Eq. (6)-(9)], we 
quoted the data of remnant mass (ui(m,Z)), metal mass 
(mz(m, Z)) and dust mass (m,d(m, Z)) of stars with mass 
m and metallicity Z from some of previous works. 

In this work, we consider AGB stars and SNe II as stellar 
sources, but neglected the SNe la for simplicity. Nozawa et 
al. (201 1) recently proved that SNe la produce little amount 
of dust. Further, Calura, Pipino & Matteucci (2008) showed 
that the dust destruction rate by SNe la is about 1/10 of that 
by SNe II. As for the metals ejected by SNe la, they play 
an important role in the chemical evolution of galaxies (e.g., 
Matteucci et al., 2009). However, since we discuss not the 
abundance ratio of each metal but the total metallicity, we 
did not take into account the contribution of SNe la. 



In this paper, we assume that the mass ranges of AGB stars 
and SNe II are 1-8 M Q and 8-40 M Q , respectively. Also, we 
assume that all stars with initial masses m > 40 M© evolve 
to black holes without SN explosions (Heger et al., 2003). 

As for the remnant and metal masses, the data are taken 
from van den Hoek & Groenewegen (1997) for AGB stars 
with mass range 1-7 M© and metallicities Z = (5.0 x 
10- 2 , 0.2, 0.4, 1.0) Z© and from Woosley & Weaver (1995) 
for SNe II with mass range 12-40 M© and metallicities 
Z = (5.0 x 10~ 2 ,0.1, 1.0) Z©. As for the dust mass, the 
data is taken from Zhukovska, Gail & Trieloff (2008) for 
AGB stars with mass range 1-7 M© and metallicities Z = 
(5.0 x 10- 2 ,0.1,0.2,0.4,0.75,1.0) Z© and from Valiante 
et al. (2009) for SNe II with mass range 12-40 M© and 
metallicities Z = (5.0 x 10~ 2 , 1.0) Z© which are quoted 
from Bianchi & Schneider (2007). 

Although Stardust yields are not completely understood, 
theoretical predictions of SNe II recently show a good 
agreement with observations of nearby supernova remnants 
(SNRs) (e.g., Nozawa et al., 2010). We considered the cur- 
rent model based on these latest results. Yet some problems 
still remain unsolved (e.g., nucleation efficiency). As for the 
dust yield of AGB stars, we adopt similar Stardust yields of 
Valiante et al. (2009, 201 1) and Gall et al. (2011a) while 
their yields may be uncertain. However, we stress that after 
the dust mass growth in the ISM becomes effective, the dust 
abundance is insensitive to the Stardust yields (Inoue, 201 1). 
Thus, although dust yields have slight uncertainties, we can 
discuss the activation mechanism for the dust mass growth 
in the ISM without ambiguity. 
2.2 Dust destruction timescale 

It is thought that SNe are the main source of the dust de- 
struction. This dust destruction process depends on various 
parameters (density and temperature of the ISM, the explo- 
sion energy of SNe, and etc.), is very complex (e.g., Jones 
et al., 1994; Jones, Tielens & Hollenbach, 1996; Nozawa et 
al., 2006). In this work, we adopt the formula presented by 
Mckee (1989). 

The timescale of dust destruction tsn is expressed as 



TSN 



M 1SM (t) 

em swept 7SN(t) : 



(12) 



where e is the efficiency of dust destruction by SN shocks, 
and is defined as the ratio of the destroyed dust to the total 
swept dust by SN shocks, m swop t is the swept ISM mass by 
a SN shock, 7sn(£) is the SN rate. In this work, we assume 
e = 0.1 (Mckee, 1989; Nozawa et al., 2006). 
The SN rate j(t) is expressed as 



7sn(*) = ( 

J rr 



40 M 



0(m)SFR(i - T m )dm. (13) 

m C ut(t)>8 M 

The range of the integration is the mass range where the 
SNe can occur (Heger et al., 2003). So, if t < T4om ? 
7Sn(*) = 0.0. 

The swept ISM mass by SN shocks m swcpt depends on 
both the density and metallicity of the ISM. In the case of 
a higher density, since the amount of materials which block 
SN blast wave is larger, the swept mass becomes smaller. 
Further, the line cooling by metals is more efficient in the 
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ISM of a higher metallicity, as a result, the swept mass 
becomes smaller. To consider these effects, we adopt the 
fitting formula derived by Yamasawa et al. (201 1) 



volume-to-surface ratio of the grains, then we obtain 



pt = 1535n 



s °' 202 l(Z/Z G ) + 0.039]" - 289 [M ] 



(14) 

where tisn is the ISM density around SNe, we assume 
«sn = 1-0 cm~ 3 as a representative value. 
2.3 Metal accretion timescale 

Dust mass in galaxies increases due to not only supply 
from stars but also accretion of refractory elements onto 
preexisting dust in clouds (e.g., Liffman & Clayton, 1989; 
Draine, 2009; Jones & Nuth, 2011). This accretion process 
is called "dust mass growth". Here, we use a term "clouds", 
so that it stands the cool component in the ISM. In our study, 
we neglect volatile dusts. Although, indeed, they exist in 
clouds, if clouds disappear or the temperature goes up, such 
dusts evaporate. More precise treatment will be shown in 
future work. The dust mass growth rate in clouds can be 
expressed as (e.g., Hirashita, 2000; Inoue, 2003, 201 1) 



dM d 



r,N7r(a 2 )apT(v), 



(15) 



where r\ is the mass fraction of the clouds, N is the number 
of dust grains, (a 2 ) is the 2nd moment of a grain size a, a 
is the mean sticking coefficient of metals, p K ^ s is the mass 
density of gaseous metals that are not contained in dust, and 
(v) is the mean velocity of metals in gas phase. Since we 
assume a spherical dust grain for simplicity, we have 



4ir(a 3 )a 



N 



M A _ 3M d 
77id At: (a 3 ) a' 



(16) 



(17) 



where md is the mean mass of the grain, (a 3 ) is the 3rd 
moment of a grain size, and a is the mass density of solid 
matter within the grains. Also, 



(18) 



where Pjf M is the averaged mass density of the interstellar 
clouds where the accretion process occurs and S is the dust 
abundance in the metal mass. Thus, defining the accretion 
timescale as 



4(a 3 )cr 



Ha 2 )ap\i M Z(v) ' 



(19) 



we obtain the fourth term in right hand side of Eq. (4). Also, 
the mass density Pjf M is estimated in terms of the hydrogen 
number density, rijj, as pfg M = pm-RUn, where p is the 
mean atomic weight (we assume p = 1.4; i.e. the ratio of the 
number of a hydrogen atom and a helium atom is 10 : 1 in 
clouds), and mn is the mass of a hydrogen atom. 

We assume that a = 1.0 (i.e. if a molecule/atom collides 
a dust grain, it sticks certainly) and a = 3 gem -3 (sili- 
cate). Considering that the dust mass growth depends on the 



r acc a 2.0 x 10 7 
a 

x 



0.1pm J V 100 cm 



»H 



T 
50 K 



0.02 



— TaccO % 



-1 



[yr] 

(20) 



where r acc ,o = 4.0 x 10 yr, and we adopted a = 0.1 pro., 
?ih = 100 cm~ 3 , and T = 50 kQ. The typical size of grains 
a is defined as (a 3 ) /(a 2 ) in Hirashita & Kuo (201 1). 

As mentioned above, we conservatively adopt a = 0.1 pm 
as a fiducial value (e.g., Inoue, 2011). Small grains may 
be depleted by coagulation in molecular clouds (Hirashita 
& Yan, 2009), which strengthens the importance of large 
grains. The importance of large grains is further enhanced 
given that the grain size distribution tends to be biased to 
large (a ~ 0.1 pm) size by the destruction within SN rem- 
nants (Nozawa et al., 2007). Thus, we assume a ~ 0.1 pm to 
estimate the dust mass growth timescale. Although we basi- 
cally adopt a — 0.1 pm, we also examine a = 0.01 pm for a 
quick growth case later. Indeed, the MRN grain size distribu- 
tion (Mathis, Rumpl & Nordsieck, 1977) has a = 0.01 pm 
(Hirashita & Kuo, 2011). In reality, the grain size distribu- 
tion in galaxies changes with time due to some processes 
(e.g., SN shocks, accretion, and etc.). As for the contribution 
of the evolution of the grain size distribution, we have been 
preparing Asano et al. (2012). 

In this paper, we adopt only rj = (no accretion growth) 
or 1 in order to avoid any fine-tuning. In fact, the effect of a 
different choice of r\ can be offset by a different choice of nn 
and T. This allows us to merge uncertainties of r\, nn, and 
T into the value of r acCi o. We set r acCj o = 4.0 x 10 5 yr as 
a fiducial value (e.g., Inoue, 2011). Other choices of r acc .o 
result in a different timing of the growth activation. This is 
explicitly expressed in Eq. (27) later. 

3. What drives dust mass growth in the ISM? 

In this section, we investigate what determines the point 
where the grain growth in the ISM becomes effective to the 
total dust mass in galaxies. 

3.1 Contribution of each physical process to the total 
dust mass in galaxies 

To examine the point where the dust mass growth becomes 
efficient as a main source of dust in galaxies, at first, we 
compare each process of dust production. In Fig. 1, we 
show the contributions of stars (solid and dotted lines), dust 
destruction (dashed line) and the dust mass growth (dot- 
dashed line) to the total dust mass in a galaxy. Solid, dotted, 
dashed and dot-dashed lines represents the contributes of the 
1st, 2nd, 3rd and 4th terms in r.h.s. of Eq. (4), respectively. 
The total baryon mass M to t (the sum of the total mass of 
stars and ISM in a galaxy) and r\ are 10 10 M Q and 1.00, 
respectively. However, as mentioned in Sect. 2.1, since M to t 
is just a scale factor, M to t just changes the values of the 
contributions of these processes linearly. 



1 This temperature corresponds to (v) = 0.14 km s _1 . We assume 
Amji(v) 2 = kT and adopt A = 20 (Amu is the mean mass of the 
colliding atoms) (Spitzer, 1978). 
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Fig. 1 . Time evolution of each dust production/destruction rate [terms in rig ht hand side of Eq. (4)] with M tot = 10 10 M and r] = 1.00. The star 
formation timescales are set to be 0.5 Gyr (top-left panel), 5 Gyr (top-right panel) and 50 Gyr (bottom panel). Solid, dotted, dashed and dot-dashed 
lines represent the rates of dust injection into stars, dust ejection by stars, dust destruction by SN shocks and dust mass growth in a galaxy, respectively. 



From these figures, we find that although the contribution 
of the ejected by stars is the biggest in the early stages, as 
the time passes, the main contributor of the dust produc- 
tion switches to the dust mass growth in the ISM at a point 
("switching point")- For example, Liang & Li (2009) pointed 
out that dust produced by SNe II is predominant to the dust 
budget in galaxies at high-z Universe (z > 5) using the ex- 
tinction curves of GRB host galaxies at high redshifts. Their 
results are in good agreement with this work. Further, the 
process of the dust mass growth is expected to explain the 
dust amount in the Milky Way or dusty QSOs at high red- 
shifts (e.g., Zhukovska, Gail & Trieloff, 2008; Draine, 2009; 
Michalowski et al., 2010a; Valiante et al„ 2011). So, what 
determines the switching point? We will discuss the point in 
next subsection (this is the main topic in this paper). 

After the dust mass growth took place, the contribution of 
dust destruction by SN shocks approaches that of the dust 
mass growth. Thus, after the dust mass growth becomes 
efficient, dust amount in galaxies determines the balance 



between the contribution of dust destruction and that of the 
dust mass growth in the ISM (see also Inoue, 201 1). 

In addition, we also observe that the increase of the contri- 
bution of the dust mass growth (dot-dashed line) has a peak, 
after that, the increase slows down. In other words, the dust 
mass growth becomes ineffective. We consider the reason of 
this. In Fig. 2, we show the time evolution of S (= Mz/M^) 
with various star formation timescales. We find that after 
the value of S increases rapidly, the value saturates. Hence, 
after the dust mass growth becomes effective, most of met- 
als form dust. Thus, the reason why the dust mass growth 
becomes ineffective is the depletion of metals. 

As shown in Fig. 2, the values of 6 for all tsfs converge 
to ~ 1. In contrast, the value for the Milky Way is about 
0.5. However, since it can be adjusted by adopting different 
r\, we do not try to fine-tune the convergence value of 5 in 
this study. Inoue (2011) showed that the convergence value 
of S is determined by the balance between the contribution 
of dust destruction by SN shocks and that of the dust mass 
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Fig. 2. Time evolution of the fraction of metals in dust, <5 with r\ = 1.00. 
Solid, dotted and dashed lines represent rgp = 0.5, 5, 50 Gyr, respec- 
tively. Cross symbols mark the switching point for each rgp . 



growth (for details, the product of r acCj o and er7i swep t). 
3.2 Critical metallicity for dust mass growth 

In this subsection, we introduce the main topic of this 
paper, the critical metallicity. This is a metallicity at the 
switching point (see Sect. 3.1). In our model, the sources 
of dust are stars (AGB stars and SNe II) and the dust mass 
growth in the ISM. 

In order to find the critical metallicity, we compare the 
second term with fourth term in right hand side of Eq. (4). 
First we consider the second term [Eq. (9)]. From Eq. (9), if 
D is defined as 



100 M 

md(m, Z(t — T m ))<p(m)dm = D, 

'nicut(t) 

then Eq. (9), with Eq. (5) can be approximated as 

, Mism 



Yd — D- 



TSF 



(21) 



(22) 



This equation is exact if SFR is constant. 

Next, we consider the fourth term in right hand side of 
Eq. (4). Since r acc = T acCj0 Z^ 1 and M d = 8M Z = 
SZAIism, the dust mass growth term is 



q- 



M d {l-S) _ r)5(l-S)Z 2 M ] 



ISM 



7"acc,0 



(23) 



From Eq. (22) and (23), the metallicity at which the in- 
creasing rate of dust mass by the dust mass growth exceeds 
the dust production rate by stars is, then, presented as fol- 
lows: 



Z 



D 



7"acc,0 
TSF 



(24) 



r)5(l-5) 

Thus, if the metallicity of a galaxy is larger than above metal- 
licity, we should consider the effect of dust mass growth in 
the galaxy. Here, we refer to the metallicity as the criti- 
cal metallicity Z CI , which is the metallicity at the switch- 
ing point. To obtain the value of Z cx , hereafter, we adopt 



8 = 0.02 and D = 5 x 10~ 4 . As for the value of 8, from 
Fig. 2, the value of 8 ranges 0.01-0.04 at the switching point 
for each tsf- Further, although 8 is dependent on time, be- 
fore the dust mass growth becomes effective to the total dust 
mass, the 8 is determined only by the contribution of stars 
(see Appendix A). This contribution of stars ranges 0.01- 
0.04 in our calculation (Fig. 8 in Appendix A). Thus, we 
adopt (5 = 0.02 as a representative value. Also, since we 
found from numerical calculation that the range of D is 10~ 4 
to 10~ 3 , we adopt D = 5x 10~ 4 as a representative value. 

For the reader's convenience, we compare our discussion 
with a similar work by Inoue (2011). Inoue (2011) defined 
a critical metallicity to compare the contribution of the dust 
destruction by SN shocks with that of the dust mass growth. 
Thus, the critical metallicity in Inoue (2011) is a metallicity 
which the contribution of the dust mass growth exceeds that 
of the dust destruction. In contrast, our critical metallicity is 
a metallicity which the dust mass growth becomes the main 
source of the increase of dust (the contribution of dust mass 
ejected by stars is the main source of dust at early stage of 
galaxy evolution). Readers who are interested in both works 
should keep this difference in mind. 

In Fig. 3, we show the relation between the metallicity 
and dust-to-gas mass ratio (left panel is normalized by the 
critical metallicity, while right panel is not normalized) for 
t"sf = 0.5, 5, 50 Gyr. In Fig. 3, we find different evolution- 
ary tracks depending on tsf in right-hand panel, whereas 
these tracks with metallicity normalized to the critical one 
are well overlaid on each other in the left-hand panel. This 
clearly shows that the critical metallicity is truly essential: 
The dust mass growth becomes efficient not when the galac- 
tic age reaches a certain value but when metallicity exceeds 
the critical value. In the left panel, the dust-to-gas mass ratio 
increases rapidly after metallicity exceeds the point of the 
critical metallicity. Thus, dust produced by the dust mass 
growth becomes efficient to the dust mass of a galaxy if 
Z > Z cr in the galaxy. 

In Fig. 4, we show the critical metallicity as a function 
of tsf- As shown in the figure, the critical metallicity be- 
comes larger if the star formation timescale is shorter. As 
mentioned in Sect. 1 and 3.1, the dust mass growth is ex- 
pected to be the dominant source of dust in various galaxies 
(e.g., the Milky Way: Draine (2009) and dusty QSOs at high 
redshift§ Valiante et al. (2011), among others), in spite of 
different star formation timescales and ages of the galaxies. 
In terms of the critical metallicity, we can explain the reason 
of it in a coherent manner; metallicity in the galaxies just 
exceeds the critical one. 

Here we emphasize the importance of the metallicity de- 
pendence in the accretion growth timescale, r acc . Pipino et 
al. (2011) argued that the dust mass growth was important 
to explain the observed huge mass of dust in high-z QSOs. 
However, they did not seem to consider the metallicity de- 
pendence in T ac( 0. Fig. 5 shows the effect of the metallic - 



2 It is considered that QSOs in high-z Universe have high star formation 
rate and even and larger subsolar metallicities (e.g., Matsuoka et al., 2009; 
Juarez et al, 2009). 

^Indeed, their equations (4) and (5) have the metallicity dependence. How- 
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Fig. 4. Critical metallicity as a function of rgp with r\ = 1.00. 



ity dependence on the dust mass evolution. When we con- 
sider the dependence properly, the timing where the dust 
mass growth becomes effective is delayed until the metal- 
licity exceeds the critical value as discussed above. On the 
other hand, if we omit the dependence and adopt a constant 
value for r acc , the timing of the growth becomes effective on 
the total dust mass is determined by just the adopted r acc . 

In the above discussion, we have focused on the critical 
metallicity. One may be, however, interested in its relation 
to the time, t CI , which is a galactic age when the metallicity 
in a galaxy reaches the critical metallicity. Here, in order to 
understand the importance of the dust mass growth in various 
galaxies with various star formation timescale more clearly, 
we demonstrate the relation between the critical metallicity 



Fig. 5. Time evolution of the total dust mass with (solid line) and without 
(other lines) metallicity dependence for the dust mass growth timescale. 
Here we adopt t sf = 5 Gyr, M to t = 10 10 M Q and 77 = 1.00. The 
three values of r aC c represent Z = 0.2 Zq (dotted line), Z = 1.0 Zq 
(dashed line) and Z = 2.0 Zq (dot-dashed line), respectively. 



Z cr and the time t CT . However, we stress that the metallicity 
is more fundamental because t CI is determined by the critical 
metallicity. 

We derive the relation with using Eq. (2) and (3). If A is 
defined as 



100 M 



mz(m, Z(t — T m ))(j)(m)d j m = A, 



'm cu t(t) 

then, the relation is expressed as 

TSF 



(25) 



(26) 



ever, their adopted timescales in table 1 seem to omit the dependence finally. where A is a constant, and is about 0.018 in our calcula- 
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Fig. 6. t cr as a function of rgp with r/ = 1.00. 

tion. If inflow process occurs, since dust-to-gas mass ratio 
and metallicity become smaller than the case without inflow, 
it is sufficient to consider the larger tsf- Figure 6 shows 
t CI as a function of tsf- From this figure, we find that t CI 
becomes shorter if tsf is shorter, while in the case of Z CI , 
the trend is opposite (Fig. 4). This is explained by following 
reason. If the t$f is short, the fast progress of the star for- 
mation makes the metallicity become large at an early stage 
of galaxy evolution. Hence, although the critical metallic- 
ity is large in the case with short tsf, t CI is short due to the 
fast evolution of metallicity. From the relation between z cr 
and t cr , we can understand that the dust mass growth is the 
main dust production of the Milky Way (evolved galaxy) and 
dusty QSOs (in high-z Universe). 

Next, we discuss the effect of the grain size distribution. 
Recently, Hirashita & Kuo (201 1) showed that the dust mass 
growth in the ISM depends on the grain size distribution of 
a galaxy. To examine the effect of the grain size on the dust 
mass growth, we consider the case with a = 0.01 /im (the 
case with a = 0.1 /im is the fiducial case in this paper). 
Figure 7 shows the relation between metallicity and dust-to- 
gas mass ratio for a = 0.1 /im (solid line) and 0.01 /im 
(dotted line) with tsf = 5 Gyr and rj = 1.00. From 
Eq. (20) and (24), the critical metallicity as a function of a is 
expressed as 

Z cr (a) = (a/0.1 /im) 1/2 Z cr (0.1 fim). (27) 

From Fig. 7, we find that the evolutionary tracks with dif- 
ferent a show the almost same behavior if we introduce the 
critical metallicity for each value of a. Thus, although Hi- 
rashita & Kuo (2011) showed that the critical metallicity is 
sensitive to the grain size distribution, the mechanism that 
the critical metallicity determines the timing on which the 
dust mass growth becomes dominant to the growth of the to- 
tal dust mass in a galaxy does not change. As for the dust 
evolution considered the evolution of the grain size distribu- 
tion (including the effects of the stellar dust, SN destruction 
and accretion) in a galaxy, This issue will be extensively dis- 
cussed in our next work (Asano et al., 2012, in preparation). 



Fig. 7. The relation between metallicity and dust-to-gas mass ratio for 
a = 0.1 /im (solid line) and 0.01 /im (dotted line) with r$p = 5 Gyr 
andr? = 1.00. FromEq. 27, Z cr (0.01 fim) ~ 0.32 Z cr (0.1 fim). 



4. Conclusions 

In this work, we constructed a galaxy evolution model 
taking into account the metallicity and age dependence on 
the various dust sources (AGB stars, SNe II and growth in 
the ISM) to investigate what is the main driver of the grain 
growth which is expected to be the dominant source of dust 
in various galaxies with various star formation timescales. 

We have found that the timing that the dust mass growth 
in the ISM becomes effective is determined by metallicity. 
If metallicity in a galaxy exceeds a certain critical value, 
critical metallicity, the dust mass growth becomes active 
and the dust mass rapidly increases until metals are depleted 
from the ISM. This critical metallicity is larger for a shorter 
star formation timescale. The dust mass growth is thought 
to be the dominant source of dust in evolved galaxies like 
the Milky Way and young but dusty and massive QSOs at 
high redshifts. The importance of the dust mass growth in 
such a diversity of galaxies can be explained clearly in terms 
of the critical metallicity; the dust mass growth in the ISM 
is regulated by metallicity, and we stress that the critical 
metallicity works as an indicator to judge whether the grain 
growth in the ISM is dominant source of dust in a galaxy, 
especially because of the strong and nonlinear dependence 
on the metallicity. 
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Appendix A. What dominates S before the dust 
mass growth becomes effective? 

Figure 8 shows the time evolution of the ratio of dust mass 
produced by stars to metal mass ejected by stars, 5 sta r (thick 
lines) and 5 (thin lines) with rgp = 0.5, 5, 50 Gyr. From 
Eq. (7)-(9), the ratio <5 s tar is expressed as 



r(*) = 



Yd® 



Rz{t) + Y z {t) 



(A.l) 



From Fig. 8, we observe that the evolutionary tracks of <5 sta r 
are the almost same tracks before the dust mass growth be- 
comes effective to the total dust mass in galaxies. Hence, 8 
before the dust mass growth becomes effective is determined 

by <5 sta r- 
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